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Abstract — This paper describes the design of a modified 
tent map characterized by a uniform probability density 
function. The use of this map is proposed as an alternative 
to the tent map and the Bernoulli shift. It is shown that 
practical circuits implementing the latter two maps may pos- 
sess parasitic stable equilibria, fact which would prevent the 
desired chaotic behavior of the system. On the other hand, 
commonly used strategies to avoid the parasitic equilibria 
onset also affect the uniformity of the probability density 
function. Conversely, the use of the proposed tailed tent 
map allows to assure a certain degree of parameter devi- 
ation robustness, without compromising on the statistical 
properties of the system. 



I. Introduction 

Discrete time chaotic systems based on simple, one- 
dimensional, piecewise linear maps can be conveniently im- 
plemented using current mode operation Low circuit 
complexity makes them appealing for many applications, 
like secure communication schemes [2], noise generators 
and analog random number generation for stochastic neu- 
ral models [sj. Design requirements can be posed either 
on the spectrum or on the probability distribution of the 
chaotic signal. For instance, the random number generator 
needed in |3] should produce samples characterized by a 
uniform probability density function (PDF). 

Although systems based on the tent map (Fig. [l]A.) and 
on the Bernoulli shift satisfy this requirement, their imple- 
mentation fl], I4] is not completely straightforward, due to 
the possible onset of parasitic stable equilibrium points, as 
result of the unavoidable errors introduced by the physical 
realization of the electronic devices. Since this fact would 
not allow the desired chaotic behavior, special strategies 
must be adopted to make the system robust to parameter 
deviation. Unfortunately, commonly adopted solutions [4] 
also prevent the achievement of a uniform PDF. 

The aim of this paper is to show that a one-dimensional 
system based on a modified version of the tent map, re- 
ferred as tailed tent map (Fig. [lj3), possesses both a uni- 
form PDF and a chaotic behavior inherently robust to im- 
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Fig. 1 

The tent map (A) and the tailed tent map (B). Normalized 

GRAPHS IN / = [0, 1]. 



plementation errors, being better suited for carrying out, 
for instance, the random number generator needed in [3]. 
In order to show the feasibility of the system, a design in a 
standard 1/im, n-well, single-poly CMOS technology is also 
proposed. Several circuit simulations confirm the uniform 
distribution of the samples produced by the system. 

A. Mathematical background 

Let / = [a, b] be an interval of the real line TZ and con- 



sider the one-dimensional dynamical system 



M{Xr. 



(1) 



where M is a nonsingular transformation of / into itself. 
When designing a pseudo random number or a noise gener- 
ator based on system ([T]) two properties are particularly im- 
portant, namely i) chaotic behavior and ii) uniform prob- 
ability density. 

i) Chaotic behavior is assured when the system ([T]) has a 
positive Lyapunov exponent, defined as [Sj: 



T-l 



lini - Vln|Af'(xi)| 



(2) 



1=0 



Roughly speaking, the Lyapunov exponent is an index of 
the sensitivity of the system to initial conditions and there- 
fore of its level of unpredictability. From ([2]) it immediately 
follows that systems based on everywhere expanding maps 
(for which infa;^/! Af'(a;)| > 1) are chaotic. 
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ii) Consider an interval 5* C / and let xs be the relative 
characteristic function. Then, the "average time" spent in 
S by an orbit originating in xq may be defined as 



lim — 

n— ^oo 71 



n-1 

E 

i=Q 



(3) 



where is the «-th iterate of the map. Note that if limit 
([3]) exists, it depends on the choice of the initial condition 
Xf). Birkhoff ergodic theorem fs' gives a sufficient condition 
for the existence of limit ([3| "almost" independently from 
xq. In particular, it states that, if A/ is measure preservinj^ 
and ergodic then there is a unique PDF p for M such that 



n-l „ 

lim - ^ xs(MXa;o)) = / p{m 



(4) 



for any measurable S* C /. Basically, ergodicity provides 
that any initial condition randomly chosen in / leads to a 
trajectory characterized by the same statistical properties. 
Furthermore, ergodicity enables extraction of information 
about time-based statistics from the a priori knowledge of 
the PDF. In particular, it can be shown that, if M is er- 
godic, the PDF p is the unique invariant under the Perron- 
Frobenius operator ^ Vm : ^ defined as 

n/wm= / piOd^ v^c/, (5) 

namely Vm [p] = P- 

In the following, we will consider the design of the tailed 
tent map in order to obtain an ergodic system and we 
will thus determine its uniform PDF as a fixed point of 
Vm- Eventually, we will verify its robustness to break- 
points misplacement both via analytic considerations and 
by approximating p using time series data extracted from 
simulations. 

II. Systems based on the tent map 
Consider a tent map based system ([T]), namely assume 



[0, 1] and 



M{x) = l-2\x- 1/2|. 



(6) 



Since the map is everywhere expanding, the system is 
chaotic. Moreover, it can be proved that the system is 
ergodic and characterized by a uniform PDF [5j. 

Fig. 2] shows a suitable circuit implementation for the 
map p Transistors Mtp, Mtn and Inv act as an active rec- 
tifier: when lin is positive (negative) Mjp (Mtn) is on and 
Mtn {Mtp) is off. Since current mirrors CM„i and CMpi 
have a 2:1 ratio, one obtains lout = Iref — 2|/i„|. The sim- 
ulated static characteristic is represented in Fig. [Sj'V and 
is the same as ^ apart from an uninfluential change of 
the reference system. Considering the implementation of 
equation ([T]) by a current mode circuit, the necessary de- 
lay operation can be realized by cascade connecting two 
switched current (SI) sample and hold stages with comple- 
mentary phase clock \V . 

^ Af is said to be measure preserving if a measure exists such that 
fi(M~'^ (A)) = fi{A) for any measurable set AC I. 
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Fig. 2 

A CIRCUIT FOR THE TENT MAP. 



A. Parasitic stable equilibrium and strategies to avoid it 

When the map is implemented using an analog circuit, 
some misplacement of the breakpoint Bp is unavoidable. 
If the ordinate of the breakpoint exceeds its assigned value 
due to a slope error, a parasitic stable equilibrium is intro- 
duced. In order to see how this happens, it is necessary 
to consider the circuit static characteristic outside the map 
domain. Due to the fact that the current mirrors or the 
current rectifier transistors (Fig. [2]) reach their maximum 
allowable current, the map slope decreases and flattens, 
as shown in Fig. [Sj'V. As a consequence a parasitic equi- 
librium point exists (point E in Fig. [s]) characterized by 
M'{xe) < 1, i.e. the equilibrium is stable. Although this 
equilibrium always exists, it ideally lies out of the map 
invariant set. On the contrary, if the breakpoint ordi- 
nate is greater than its assigned value (independently of 
the magnitude of the misplacement), the map invariant set 
abruptly expands, including also point E. Therefore, af- 
ter a transient phase, the system trajectory will converge 
towards it, as schematically represented in Fig. |3j3. The 
abrupt expansion of the invariant set is due to the fact that 
its closure contains the unstable equilibrium point F. As 
soon as the invariant set grows to include F in its interior 
it must also contain the whole interval [xE,xp]. 
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I invariant set when brk.point is too high , 

(B) 
Fig. 3 

(A): Existence of a parasitic stable equilibrium and (B): 

CONVERGENCE OF THE SYSTEM TRAJECTORY TOWARD IT, DUE TO A 
MISPLACEMENT OF Bp IN THE MAP. 



Note that besides breakpoint misplacement, also noise 
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added to the system state may allow the system trajectory 
to converge to E. 

In order to avoid this problem, several strategies can 
be used. All of them are based on small modifications of 
the majj^ and are equivalent to designing a system whose 
breakpoint ordinate is less than 1, to assure an adequate 
margin. Two typical solutions are shown in Fig. [4j 





(A) 



Fig. 4 

Two STRATEGIES TO AVOID CONVERGENCE TO THE PARASITIC 
EQUILIBRIUM. (A) SHOWS A STRATEGY WHICH AFFECTS THE SLOPE OF 
BOTH SEGMENTS OF THE TENT MAP, WHILE (B) SHOWS A STRATEGY 

WHICH AFFECTS THE SLOPE OF THE SECOND SEGMENT ONLY. By 
RESCALING THE GRAPH IT CAN BE SEEN THAT ALSO THE CASE (B) IS 
EQUIVALENT TO LOWERING THE BREAKPOINT. 



The margin one must allow mainly depends on three 
factors: 

1. the maximum random error one can expect for the cir- 
cuit static characteristic. This originates mainly from tran- 
sistors mismatching; 

2. the error introduced by the two SI sample and hold per- 
forming the delay operation. This is mainly due to clock 
feedthrough and appears as an undesired signal added to 
the system state; 

3. the error due to the dynamic response of the current 
comparator used to implement the breakpoint (Fig. [2]). 
This error manifests itself as a sort of hysteresis and ap- 
pears when the circuit has to switch immediately over or 
below threshold. 

In particular, if speed is an important requirement, tran- 
sistors dimensions must be reduced in order to contain par- 
asitic capacitances effects and the first factor will have great 
importance. On the other hand, if a low-power circuit is 
desired, one has to accept a worsening of the current com- 
parator dynamic response and a consequent increased rel- 
evance of the third factor. Therefore, particularly in fast, 
low power systems the required margin can be appreciable. 

B. Pitfalls in the strategies to prevent convergence to par- 
asitic equilibria 

Adopting one of the above mentioned strategies implies 
sacrificing PDF uniformity. For instance. Fig. [5] shows the 
consequences on the PDF of modifying the map as shown 
in Fig. |4]A_. 

^ Without loss of generality we will refer our considerations to nor- 
malized maps defined in / = [0, 1]. 
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Fig. 5 



Consequences on the PDF of the strategy shown in Fig. |4]4. 



The graphs are obtained simulating system ([T]) for a tent 
map in which the ordinate of Bp has been decreased by 2%, 
i.e. for a system tolerant to breakpoint misplacements up 
to 2%. In Fig. [5]A. it is represented the case in which no 
misplacement occurs (ideal case), while Fig. [5^3 shows the 
worst case in which a 2% misplacement takes the break- 
point even further down. 

By similar considerations, it can be seen that systems 
based on the Bernoulli shift suffer from the same inconve- 
niences mentioned above 

III. The new approach: systems based on the 

TAILED TENT MAP 

The tailed tent map (Fig.[lj3) is a parametric map, whose 
secondary breakpoint also defines the position of the main 
one, namely 

{l-Xb2 + 2x X e [0, a::fc2/2) 
1+Xb2-2X X e [xb2/2,Xb2) , (7) 
1 - X X e [Xb2, 1] 

with 1/2 < Xb2 < 1. In the following, it will be shown that 
system ([I]) with / = [0, 1] and M given by ([t]) has the same 
basic properties as the one based on the tent map, yet is 
more robust to implementation errors, whenever a uniform 
PDF is the desired goal. 

A. Basic properties 

A piecewise linear map M :/—>■/ is also a Markov map 
if it takes partition points into partition points, namely a 
partition W = {(oq, ai), . . . , (ajv-i, oat)}, Uq = a, = b, 
exists such that if Q = {ao, . . . , a^} then M{Q) C Q. It is 
easy to show that whenever Xb2 is a rational number, ([t]) is 
a Markov map. Moreover, in this case, M satisfies also the 
conditions 1 — 4 of Definition 3 in [t] which assure the PDF 
unicity and the ergodicity of the map. Thus, it is easy to 
prove that 

P{x) = X[o.i\{x) (8) 

is the invariant PDF. In fact, it can be easily seen that as 
any x € I has a finite number of counterimages, then Vm 
can be also written as 



E 

M(y) = 



pjy) 

|M'(y)|- 



(9) 
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Thus, by substituting (|8| in ([9| one immediately gets 

^m[X[o,i]](2^) = 1 • X[o,i-x,2]i^) + Q + ^) • 

Knowing the invariant PDF, the Lyapunov exponent for 
the map can be computed as [s] 

\n\M'{x)\p{x)dx, (10) 

which leads to h = xi,2 In 2, proving system chaoticity. 

B. Inherent robustness of the map 

First of all, note that one can design the map in such 
a way that point E does not exist, as shown in Fig. |6]\ 
(dot-dashed line). Furthermore, even if E exists, the al- 
teration in the invariant set consequent to breakpoint mis- 
placement is now graceful, i.e. the smaller the misplacement 
the smaller the invariant set change, as long as this set does 
not include the unstable equilibrium point F. Therefore, 
up to a certain misplacement the invariant set can grow 
without including E (Fig. |6j3). More precisely, the lower 
Xb2, the higher the system tolerance. As a consequence, 
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Graceful alteration of the invariant set for the tailed 

TENT MAP. 



a system based on the tailed tent map can be designed 
without having to compromise on the PDF in principle. 

The effect of implementation errors is shown in Fig. [Tj 
In particular. Fig. [7j^ is the ideal case (i.e. the map is 
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(A) (B) 
Fig. 7 

PDF DEGRADATION DUE TO BREAKPOINT MISPLACEMENT FOR THE 
TAILED TENT MAP. 



implemented correctly) and Fig. [7^3 represents the case of 
a 2% Bpi misplacement. The superiority of the case in 
Fig. [TJ'V over the one in Fig. [sj^ is evident and also the case 
in Fig. [7f3 is more similar to an ideal uniform PDF than the 
one in Fig. [5}3, since the latter even presents a much larger 
range of values characterized by null probability density. 
The performance improvement becomes even more visible 
when one considers probability distributions rather than 
densities. 

C. Circuit CMOS implementation and performances 

Fig. [S] shows a schematic of the proposed circuit imple- 
mentation for the tailed tent map. Its left side is the same 
as in Fig. [2j while it right side allows to obtain the slope 
variation necessary to carry out the "tail." Assuming that 
current mirrors CMp2 and CM„3 have a 1 : 1 ratio, the 
output current is lout = Iref - 2|/i„| + rect(/i„ - /re/2), 
where rect(a;) = l/2(a; -I- |a::|), namely a tailed tent map 
apart from the usual change of system reference. 

Note that the extra silicon area requested to generate 
the tailed tent map is only a small part of the area used 
for the complete circuit. In fact, this includes two sample 
and hold circuits which use a large proportion of the total 
area. 





Fig. 8 

A CIRCUIT FOR THE TAILED TENT MAP. 



Fig. [9] shows the results of device level simulations 
for the circuit in Fig. [8] designed using a standard Ifim 
CMOS technology. To obtain high linearity for the map 
branches, high swing cascode mirrors have been used while, 
to improve maximum frequency operation, a careful de- 
sign of the sample and hold stage has also been consid- 
ered. Fig.[9]A_ shows a typical chaotic waveform for /o«t(fc). 
Fig. [9J3 represents the approximation of the PDF obtained 
by a post-processing executed on the transient analysis out- 
put, for a 750 KHz operating frequency. Deviation from the 
ideal line is not only due to circuit misbehavior, but also 
to the finite number of samples used to extract the PDF. 

In conclusion, the new proposed solution shows several 
improvements compared to previously published results, 
with respect to the achievable PDF [l] and parameter de- 
viation robustness [1] . 
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(A) (B) 
Fig. 9 

A SAMPLE SEQUENCE (20 SAMPLES) AND THE PDF OBTAINED 
THROUGH HSPICE SIMULATIONS AND POST-PROCESSING (50K 
SAMPLES USED FOR POST-PROCESSING). 
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